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I comment on a recent paper by Ruiz and Tsallis [Phys. Lett. A 376, 2451 (2012)] claiming to 
have found a "g-exponential" generalization of the large deviation principle for strongly correlated 
random variables. I show that the basic scaling results that they find numerically can be reproduced 
with a simple example involving independent random variables, and are not specifically related to 
the g-exponential function. In fact, identical scaling results can obtained with any other power-law 
deformations of the exponential. Thus their results do not conclusively support their claim of a 
g-exponential generalization of the large deviation principle. 



Ruiz and Tsallis [T] have published a recent a paper 
in which they claim to have found a generalized large 
deviation principle (LDP) for a process involving strongly 
correlated random variables. This LDP supposedly takes 
the form not of an exponential, as in standard large devi- 
ation theory [2] , but of a so-called g-exponential function, 
which is one possible generalization of the exponential 
function. My purpose here is to show that the main re- 
sults of PP, which take the form of scaling plots, are not 
specifically related to g-exponentials or correlated ran- 
dom variables, and do not conclusively point, as claimed, 
to a q-cxponential generalization of the LDP and large 
deviation theory. I illustrate these points with a simple 
example that shows how one might be misled into defining 
"g-exponential LDPs" for distributions that do not have 
in fact the form of a g-exponential. 

The example is a classic one in large deviation theory 
showing that not all sums of random variables have an 
LDP [2, 3J. Consider the sample mean 
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and assume that the -XVs are independent random vari- 
ables identically distributed according to the Cauchy dis- 
tribution: 
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where 7 is a positive parameter that controls the width of 
the distribution. At this point I should say that, since the 
Xf8 are assumed to be mutually independent, one might 
be tempted to discard this example as being unrelated 
to the results of Ruiz and Tsallis pQ. However, it will 
become clear that the properties of S n are in fact relevant 
to discuss their claims. 

Since Cauchy random variables are stable [J] , the prob- 
ability density of S n is also a Cauchy distribution: 
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With this result we can see that p(S n — s) concentrates to 
the Dirac-delta function S(s) in the limit n — > 00 whenever 



/3 > 1. Hence the density of S n concentrates in some limit, 
which is the first ingredient that we need to define an 
LDP. However, since 
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as n — > 00, the concentration to the delta function is not 
exponential in n, which implies that p(S n — s) does not 
satisfy an LDP. In fact, if we try to define an LDP for S n 
using the standard large deviation limit, 



I(s) = lim \np(S n = s) 
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we find that the rate function J(s) is for all s. This 
means that S n does not satisfy an LDP: one could say that 
it has a zero rate function, but this gives us no information 
about S n apart from knowing that the concentration of 
p(S n ) is slower than exponential. The same conclusion 
applies to any sums of (symmetric) Levy-stable random 
variables. 

From this it is tempting to change the definition of the 
LDP in the following way. Since we know that p(S n — s) 
decays to zero for s / as a polynomial in n, we could 
try to quantify this decay by a rate function defined on 
a scale different from logarithm. This is in essence the 
proposal of Ruiz and Tsallis [T]. In particular, we could 
replace, as they do, the normal log with the so-called 
q-logarithm function defined by 



x x ~ q — 1 
ln q (x) = ' — _ - , 

which converges to the normal logarithm for q — > 1. 
For the Cauchy example, it is easy to see that 
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In, p(S n = s) - 
as n — > 00, so we can choose 
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to obtain 
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FIG. 1: Plot of the g-log of the probability density p{s) of S n 
for P = 1.5 as a function of n. Different lines correspond to 
different values of s (from to 2 in steps of 0.1, from top to 
bottom). The plot is for q = 3, the only value of q for which 
this plot yields straight lines. 
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FIG. 2: Plot of the «-log of the probability density p(s) of S n 
for ft = 1.5 as a function of n for different values of s (from 
to 2 in steps of 0.1, from top to bottom). The plot is for 
k = 2, the only value of k for which this plot yields straight 
lines. 



where c s is some constant that depends only on s. This 
behavior is illustrated in Fig. [I] with j3 = 1.5 for which 
g = 3. Any value of q not given by Eq. ^ will yield a 
plot of hi q p(S n — s) that is not linear in n, as can easily 
be checked analytically and numerically. 

If we are to follow Ruiz and Tsallis pQ , then the straight 
lines seen on the g-log plot of p(S n — s) should imply 
that this density has the form of a ^-exponential, 

e% = [l + (l-q)x]^, (10) 

the inverse of the g-logarithm. But, as can be seen from 
Eq. ((3), this is an incorrect conclusion: p(S n = s) does 
not have the form of a ^-exponential, although it has the 
same tails as this function, so the straight lines that we 
see are not a conclusive sign of a g-exponential LDP 
they are only an indication of power-law tail behavior. 
Any density that has power-law tails will lead to straight 
lines over a wide range of values when plotted on a g-log 
scale, and this even if its center does not behave as a 
power-law, as is the case for the Cauchy sum. Moreover, 
it will show straight lines when plotted with any other 
power-law deformations of the exponential, so there is 
nothing special about the q-log and the g-exp. 

To illustrate this point, consider the so-called k- 
logarithm introduced by Kaniadakis [5], 

ln {«} a; = Yk ' (U) 

which recovers the normal log as k — > 0. This deformation 
of the log is obviously different from the q-\og; yet it can 
be used in the same way as before to transform p(S n = s) 
into a seemingly linear function of n. In this case, we 
have to choose K = q — 1 to obtain 

ln{ K } p{S„ = s) ~ -d s n (12) 

as n — > oo. This is shown in Fig. [2] 

From this it is clear that if Ruiz and Tsallis [T] had 
decided to use the K-log instead of the g-log, they would 



have obtained similar results, i.e., straight lines that are 
neither necessarily related to the q-exp nor the K-exp, the 
inverse of the K-log. Therefore, their claim that they have 
found a generalization of the LDP linked specifically to 
the g-exp is misleading: the only conclusion that follows 
from their plots is that they study a distribution with 
power-law tails, having the same asymptotic form as any 
power-law deformation of the exponential. 

Of course, the reason why the K-log does the same job 
as the q-log is that they both behave as inverse power-laws 
for small arguments, so that a function with power-law 
tails must appear linear when transformed with either 
function. From this simple reasoning, any other power-law 
deformation of the logarithm (there is an infinite number 
of them) will do the job. In fact, we do not even need 
to consider deformations of the logarithm: by inserting 
p{S n = s) into 

fr(x) = y, (13) 

we can choose r = 1/(1 — j3) to obtain 

f r (p(S n = a)) ~ -n, (14) 

which is just a different way of expressing the power-law 
behavior of p(S n — s). The same trick can be applied to 
the distribution considered in [T] with the same result. 

In the end, one could try to distinguish whether the dis- 
tribution of a random variable such as S n exactly follows 
a g-exp, a K-exp or another power-law deformation of the 
exponential by studying its behavior for small values of n. 
However, this would defeat the purpose of large deviation 
theory, which is to derive general asymptotic results by 
considering specific scaling limits of stochastic processes, 
such as the large-n limit studied here, or the long-time or 
low- noise limits [2j. It is by considering these limits that 
general asymptotics are seen to arise in a broad class of 
stochastic processes. 

The results of pQ can be taken as a sign of an interesting 
power-law scaling in the tails of a probability distribution. 
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But, as I have shown, they are insufficient to conclude 
that this scaling has the form of a g-exponential or that 
it is specifically related to correlated random variables, 
as implied in pQ. Moreover, as these results only concern 
one specific process, they are obviously not enough to 
conclude that one has generalized large deviation theory 
to strongly correlated random variables. As of now, no 
such generalization exists. 
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